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Abstract. Integrated X-ray spectra of an evolving population of starburst galaxies (SBCs) are determined based on the ob¬ 
served spectra of local SBGs. In addition to emission from hot gas and binary systems, our model SBG spectrum includes a 
nonthermal component from Compton scattering of relativistic electrons by the intense ambient far-IR and the (steeply evolv¬ 
ing) CMB radiation fields. We use these integrated spectra to calculate the levels of contribution of SBGs to the cosmic X-ray 
background assuming that their density evolves as (1 -|- z)'^ up to a maximal redshift of 5. We find that at energies 10 keV 
this contribution is at a level of few percent for q < 3, and in the range of 5% — 15% for q ~ 4.5. The Compton component 
is predicted to be the main SBG emission at high energies, and its relative contribution gets progressively higher for increasing 
redshift. 

Key words. Galaxies: X-ray - Galaxies: spiral - Galaxies: star formation - diffuse radiation 


1. Introduction 

In starburst galaxies (SBGs) enhanced star formation activity 
(lasting typically ^ 10® yr) drives a chain of coupled stellar 
and interstellar phenomena that are manifested by intense far- 
infrared (FIR) emission. The SBG spectrum shows a distinctive 
large bump of thermal FIR dust emission which is correlated 
with optical emission (e.g., Silva et al. 1998), a result of the 
fact that dust emission is reprocessed starlight from a popula¬ 
tion of hot OB stars. Such stars are short-lived (10®“^ years) 
and end up as supernovae. The subsequent radio emission from 
synchrotron radiation in supernova remnants (SNRs) is tightly 
correlated with the IR emission (e.g., Condon 1992). Interest 
in SBGs stems also from the realization that these resemble 
young galaxies in the earlier universe. Indeed, a starburst phase 
was very common then, as a result of both astrophysical pro¬ 
cesses (baryonic infall and early star formation) and dynamical 
processes (close encounters and mergers which, in turn, trig¬ 
ger star formation). Consequently, the cosmic star formation 
rate (SFR) was substantially higher at epochs corresponding to 

1, with the SFR having either a peak at 2 (Madau et 

al. 1996), or a plateau out to z ~ 4 (Thompson et al. 2001). 

A primary manifestation of the starburst activity is X-ray 
emission. Given the greatly enhanced SFR, energetic phenom¬ 
ena related to stellar evolution - OB stars. X-ray binaries, 
SNRs, galactic winds, and Compton scattering of ambient FIR 
and CMB photons off relativistic electrons that are accelerated 
by SN shocks - clearly suggest that SBGs are typically more 
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powerful X-ray emitters than normal galaxies of comparable 
stellar mass (e.g., Rephaeli et al. 1991; Schmitt et al. 1997). 
The mean X-ray spectrum of SBGs reflects the diverse nature 
of high energy activity in these galaxies (Persic & Rephaeli 
2002; Rephaeli et al. 1991, 1995). 

The link between SBGs and the cosmic X-ray background 
(CXB) was made early on (Bookbinder et al. 1980; Stewart et 
al. 1982; Weedman 1987; Griffiths & Padovani 1990; Rephaeli 
et al. 1991, 1995; David et al. 1982; Ricker & Meszaros 1993; 
Moran et al. 1999). Most previous estimates of the SBG contri¬ 
bution to the CXB were essentially based on the SBG emission 
at 2 keV, with either no source evolution (Weedman 1987), or 
with some assumed degree of evolution (Griffiths & Padovani 
1990; Ricker & Meszaros 1993). Clearly, the use of the 2 
keV emission as baseline is unjustified given the wide spectral 
range (3-100 keV) over which the CXB seems to have been 
well determined (e.g., Gruber et al. 1999). The SBG contri¬ 
bution was also evaluated based on the emission over a wider 
spectral range, either in the context of an assumed spectral 
model (Ricker & Meszaros 1993), or from a statistically de¬ 
duced mean spectrum (determined from the HEAO-1 A2h-A 4 
datasets: Rephaeli et al. 1991, 1995). 

In this paper we compute the contribution of SBGs to the 
CXB incorporating recent results in the study of SBGs. In 
particular, we use the knowledge of SBG spectral properties 
gained from direct ASCA, BeppoSAX, and RXTE observa¬ 
tions of local galaxies. This, complemented with an assumed 
form for the cosmic evolution of the SBG X-ray luminosity 
density, enables us to compute the component of the CXB spec¬ 
trum which originates in SBGs. We first (section 2) review the 
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spectral properties of (local) SBGs, and proceed to discuss the 
local SBG luminosity density and its cosmic evolution in sec¬ 
tions 3 and 4. The predicted CXB spectrum is then calculated 
in section 5, followed by a discussion (section 6), and a sum¬ 
mary (section 7). Our calculations are carried out in the context 
of the Einstein-de Sitter (0^ = 1, = 0; hereafter EdS) and 

(currently favoured) flat A (0^ = 0.3, Oa = 0.7) cosmologi¬ 
cal models, with Hq = 50 km s“^ Mpc“^. 

2. X-ray spectral properties of SBGs 

Local SBGs have been repeatedly observed in X-rays (see 
Eabbiano 1989; Rephaeli et al. 1995; Dahlem et al. 1998; Persic 
& Rephaeli 2002; and references therein). In most of these 
(~ 10) galaxies the 0.5-10 keV spectra can be well fit by a 
combination of thermal, low temperature kT ~ 0.8 keV emis¬ 
sion, and photoelectrically self-absorbed emission that can be 
represented by an exponentially cutoff power law (CPE) of the 
form Tg x/kT^ 

where x is the energy in the source frame, 
r is the photon index and kT is the cutoff energy ' (Persic & 
Rephaeli 2002, Persic et al. 2002; see Eig.l). 

The low-temperature and CPE components can be inter¬ 
preted as due to emission from galactic winds and X-ray bi¬ 
naries, respectively. Erom a detailed discussion of a synthetic 
X-ray spectrum of SBGs, based on an evolutionary model of 
galactic stellar populations and on the X-ray spectra of the rel¬ 
evant emission processes, Persic & Rephaeli (2002) have sug¬ 
gested that X-ray binaries contribute most of the 2-15 keV 
emission. Both types of X-ray binaries, high- and low-mass 
systems, have spectra that can be described as variously cut¬ 
off power laws (White et al. 1983; Christian & Swank 1997). 
In general, the population-averaged spectrum emitted by a re¬ 
alistic mix of X-ray binaries (high- and low-mass systems of 
various luminosities in Galactic proportions) can be described 
as a CPE with E ~ 1.2 and kT ~ 8 keV (Persic & Rephaeli 
2002). 

The enhanced SER activity in SBGs results in a high SN 
rate and consequently - among other effects - in a higher rate 
of electron acceleration to relativistic energies. This prospect, 
coupled with the higher energy density in the enhanced EIR 
radiation field, almost certainly implies that Compton scatter¬ 
ing of the electrons by both the EIR and CMB fields yields a 
substantially higher level of nonthermal X-ray emission than 
in normal galaxies (Schaaf et al. 1989; Rephaeli et al. 1991). A 
quantitative assessment of this process in NGC253 was given 
by Goldshmidt & Rephaeli (1995). Substantial PL X-ray emis¬ 
sion may also be produced in the galactic nucleus, especially if 
the SB-driven turbulence of the gas increases the mass inflow 

* Earlier spectral studies, in which simple thermal and PL models 
were fitted to the main emission component, concluded that ^ 5 keV 
thermal or L ~ 2 PL fits were similarly successful for all analyzed 
objects (NGC253, M82: Ptak et al. 1997, Cappi et al. 1999; M83: 
Okadaetal. 1997; NGC1569: Della Ceca et al. 1996; NGC2146: Della 
Ceca et al. 1999; NGC2903: Mizuno et al. 1998; NGC3256: Moran et 
al. 1999; NGC3310, NGC3690: Zezas et al. 1998), and hence they 
were unable to reach definitive conclusions on the nature of the main 
component. 


rate onto a central black hole (e.g., Veilleux 2001). Such spec¬ 
tral components can be described as ex x^~°‘ with the photon 
index in the range a ~ 1.6 — 1.8 over a broad energy range 
(e.g., Goldshmidt & Rephaeli 1995; Rothschild et al. 1983). 

The spectral similarity shown by the local SBGs (see also 
Roberts et al. 2001) provides a reasonably good basis for a 
first approximation at modelling the population as a whole. 
We therefore assume that a typical starburst spectrum con¬ 
sists of a 0.8 keV thermal component [i.e., a bremsstrahlung 
spectrum, with the Gaunt factor numerically calculated as pre¬ 
scribed by Itoh et al. (2000)], plus a CPE with photon spectral 
index E = 1.2 and cutoff energy kT = 7.5 keV. This second 
component is photoelectrically absorbed in situ through a HI 
gas of column density A^h- The two components are normal¬ 
ized such that the thermal-to-CPL energy flux ratio in the 2-10 
keV band is 0.03. 

In addition to the wind and binary emission, we include 
also two emission components originating from Compton scat¬ 
tering of relativistic electrons by the EIR and CMB radiation 
fields. While the full calculation of the respective fluxes re¬ 
quires a detailed, self-consistent solution of the kinetic equa¬ 
tion for the electrons (taking into account their varius energy 
loss mechanisms and their propagation mode in the disk and 
halo - see, e.g., Rephaeli 1979), we greatly simplify the treat¬ 
ment by keeping only the most salient features pertinent to our 
discussion here. Assuming the typical value of 0.8 for the ra¬ 
dio (spectral energy) index of the radio synchrotron flux in the 
disk, and estimating typical energy loss times of the electrons 
by synchrotron emission and Compton scattering off the EIR 
and CMB radiation fields, we can roughly represent the latter 
two Compton components as follows: A primarily inner-disk 
EIR-scattered component with a (photon) flux that is oc 
and a (primarily) outer-disk and halo CMB-scattered compo¬ 
nent oc x~^'^ due to electrons from the higher-energy region of 
the (electron) spectrum that is steepened by radiative losses. 

The relative contribution of the Compton fluxes is substan¬ 
tially uncertain since this emission has not yet been unequiv¬ 
ocally detected in SBGs (see, e.g., Goldshmidt & Rephaeli 
1995 and Rephaeli & Gruber 2002). However, the expectation 
that Compton scattering of the enhanced relativistic electron 
population by the intense ambient EIR radiation field provides 
strong motivation for inclusion of this emission in our mod¬ 
elling of SBG X-ray emission. Our modelling of the Compton 
fluxes is based on the radio properties of the nearby SBGs 
M82 and NGC253, and well known Compton-synchrotron re¬ 
lations (e.g.. Tucker 1975) between the radio luminosity and 
deduced mean magnetic field in the disk (typically, ~ 10 /rG, 
but with appreciable uncertainty) and the predicted X-ray flux 
from Compton scattering of the radio producing electrons by 
the EIR and CMB fields. The relative strength of these two ra¬ 
diation fields largely determines (under typical conditions in 
local SBGs) the normalization between their respective fluxes. 
The energy density of the EIR radiation field, C/gr, was esti¬ 
mated by Goldshmidt & Rephaeli (1995) in the context of a 
quantitative model for the emission from warm dust. We adopt 
their result, C/gr ~ 10 Uq where Uq is the energy density in 
the CMB, as roughly typical in the inner region of a SBG disk. 
Doing so, and using the usual Compton-synchrotron relations. 
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we obtain that in local SBGs we might expect that the FIR and 
CMB Compton contributions to the overall 2-10 keV (energy) 
flux are roughly at a level of 5% and 0.5%, respectively. We do 
not yet have unequivocal evidence for nonthermal Compton X- 
ray emission in SBGs. Power-law fits to both ASCA and RXTE 
data on M82 & NGC253 were found to be formally acceptable 
(e.g.: Matsumoto & Tsuru 1999; Rephaeli & Gruber 2002), but 
a closer assessment of the residuals of such fits over the wider 
energy range afforded by RXTE leads only to what are essen¬ 
tially upper limits to nonthermal contribution. With the higher 
spatial resolution of the Chandra satellite it was clearly estab¬ 
lished that some of the 2-10 keV emission in M82 is diffuse, 
emanating from a central elliptical region, 0.8 kpc x 0.8 kpc, 
that is coplanar and coaligned with the galactic disk (Griffith et 
al. 2000). According to these authors the spectral fits to these 
measurements make a Compton origin of most of the emission 
quite likely. The luminosity of the extended emission region is 
L 2-10 — 2.2 X 10^® erg s“^ (for an assumed distance D = 3.6 
Mpc); this value constitutes ~ 5% of the total 2-10 keV lumi¬ 
nosity (T 2-10 — 4.5 X 10'*° erg s“*) measured by BeppoSAX 
(Persic & Rephaeli 2002). Given that some or most of this dif¬ 
fuse hard component in M82 is thermal as suggested by the 
presence of a substantial 6.7 keV Ee-K emission, the figure of 
5% is actually a strong upper limit on the Compton contribution 
in M82. 2 

Einally, let f{x) be the spectral profile (in the source 
frame), normalized in the band relative to which the surface 
brightness is being computed. In accord with our discussion 
above, the X-ray spectral profile of a SBG is a superposition of 
components that can be written as 

fix) = PfTBix) + -f 

+ + mfc^^ix,z) ( 1 ) 

where a: = e(l-|-z)is the energy (in keV) in the source frame 
(and x is the redshift; e is the energy in the observer’s frame), 
CTpa is the cross section for photoelectric absorption per unit 
hydrogen column density TVh within the source [in the calcu¬ 
lations we have adopted the analytic expression for apa given 
by Morrison & McCammon (1983); is taken at a nominal 
level of 1022 cm~^], and /?, pi, and 772 are constants of nor¬ 
malization that are estimated as described above. Note that the 
CMB scattered flux is z-dependent, an important fact that will 
be elaborated upon in Section 5. 

3. Local luminosity density 


NGC 253 




Fig. 1. The best-fit thermal - 1 - CPL models for M82 (left) and NGC253 
(right), based on BeppoSAX data (see Persic & Rephaeli 2002). 


Template SBG Spectrum 



Energy »>■ 

Fig. 2. The template photon SBG spectrum assumed in the present 
calculations. It consists of: (i) an unabsorbed 0.8 keV thermal 
bremsstrahlung component, (ii) an exponentially cutoff PL with pho¬ 
ton index F = 1.2 and cutoff energy kT = 7.5 keV, absorbed through 
Nil — 10^^ cm~^, and (in) a similarly absorbed PL with photon in¬ 
dex ai = 1.8. (The much fainter unabsorbed PL with photon index 
02 = 2.3 is not shown here.) The flux normalizations are as specified 
in the text. 


Sources in the IRAS Point Source Catalog with fluxes mea¬ 
sured in all four IRAS bands and with individual, identifiable 
galaxy types inhabit distinct and well-defined areas in color- 
color diagrams (Rowan-Robinson & Crawford 1989). Each one 
of these distinct components of the IRAS point-source popula¬ 
tion largely produces most of the emission at specific wave¬ 
lengths, and is characterized by a specific set of parameters for 
eq.(2o). Eor the warm-source component that peaks at ~60/im 
which is thought to be attributable to a dusty SBG population, 
the parameters in eq.(2a) take the values: 


{ C = 3.25 X 10-*Mpc“° 
logL* = 9.99 
a = 1.27 
cr = 0.626 


(25) 


In order to calculate the contribution of SBGs to the CXB we 
need to model the population at all z, beginning - of course - 
with the local source density which can be deduced from the 
EIR luminosity function. This has been determined from the 
60/im flux-limited IRAS samples (Saunders et al. 1990), 

(j)iL)d\ogL = 


^ Also for NGC 3256 (the brightest among local SBGs), Compton 
scattering has been argued to contribute appreciably to the 2-10 keV 
flux (Moran et al. 1999). 


for luminosities in the range 6 < logiL^Q/L q) < 12.5. 
Pearson & Rowan-Robinson (1996) modeled the (sub-mm, IR, 
optical) number-flux relation for extragalactic sources by in¬ 
cluding a SBG population described by a local luminosity func¬ 
tion as in eq.(2). The local 6Qpm luminosity density implied by 
eq.(2) is: 

= 8.4 X 10^ Lq Mpc“° . (3) 

Eor a sample of local SBGs observed with ASCA, the ratio of 
2-10 keV flux to 60/im flux ^ is (4 ± 1) x 10“*. Assuming this 

* The objects are: M82, NGC253, NGC2903, NGC3690. The 
60/tm fluxes are calculated using the prescription T’(60) = 2.58 x 
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to be the mean value for the whole SBG population, we have 
for the 2-10 keV luminosity density 

A-iokevW = (1.3 ±0.3) X 1037 ergs-7Mpc-3. (4) 

Note that in our previous estimates of the SBG contribu¬ 
tion to the CXB (Rephaeli et al. 1991, 1995) we did consider 
only the high end (±fir ± lO^^ Lq) of the FIR luminosity 
function as representative of the SBG population. Adopting the 
current more inclusive luminosity function (Pearson & Rowan- 
Robinson 1996) allows a more realistic estimate that accounts 
also for the emission from fainter SBGs. By integrating the 
FIR luminosity function over the above (logarithmic) range 
([6, 12.5]), we compute a density of ~ 5.9 x 10“3 Mpc“3 
for galaxies that are identified here as members of the broader 
SBG population. (Obviously, this value of the density is much 
higher than that of the more luminous galaxies, but this differ¬ 
ence is largely compensated for by the much lower mean FIR 
luminosity, Lfir ~ 10® i©, of the SBGs that we include in 
this work.) 

4. Cosmic evolution 

Galactic evolution may generally be reflected in density and lu¬ 
minosity variations over cosmic time due to the continued for¬ 
mation, interaction and merging, and through intrinsic changes 
of the luminosity due to dynamical and thermal evolution of 
stars and interstellar gas. The evolution of SBGs is likely to 
be particularly pronounced as the SB phase is thought to be 
triggered by galactic interactions and is relatively short lived. 
These imply that SBGs were more abundant at earlier times 
than would be expected based solely on pure density evolution 
of sources in an expanding background (ex (1 ± z)^). 

Observationally, the evolution of SBGs is inferred from ra¬ 
dio, IR, and optical measurements. At radio frequencies, two 
distinct populations are identified at 1.4 GHz (Benn et al. 
1993): radio-loud, and radio-faint, the first consisting of giant 
ellipticals and QSOs, and the second mainly of spirals (SBGs 
and radio-quiet QSOs). This latter faint population is respon¬ 
sible for the low-flux upturn in the differential counts, and is 
indistinguishable from the SBG-dominated 60/im population 
(e.g., Danese et al. 1987; Franceschini et al. 1988), which can 
be fit by a model using pure luminosity evolution of the form 
L{z) = L{0) (1 ± z)'^ with q = 3.1 (Rowan-Robinson et al. 
1993). At IR wavelengths, Lonsdale et al. (1990) concluded 
that the dominant population in the faint 60/im counts from the 
IRAS Faint Source Survey consists of strongly evolving SBGs. 
^ Thus, it appears that from 60/im to radio wavelengths what 
is mainly sampled is starburst activity. SBG evolution, either 

10-14j^(6o), where fi, is the nominal flux densities listed in the 
IRAS catalog (Lonsdale et al. 1985; see also Shapley et al. 2001). The 
2-10 keV fluxes are from Ptak et al. 1997 (M82, NGC253), Mizuno et 
al. 1998 (NGC2903), Zezas et al. 1998 (NGC3690). 

* In ultra-luminous IR galaxies (ULIRGs), which are the high- 
Lqoiiui tail of the ’warm’ 60/im-selected IRAS population, the star- 
formation activity (Kim & Sanders 1998a) is very strong; these un¬ 
dergo cosmic evolution (Kim & Sanders 1998b) similar to the one 
derived for the whole 60/tm-selected SBG sample. In addition, a 
preliminary XMM-Newton study of a sample of ULIRG spectra 




Fig. 3. The spectrum of the integrated background light arising from 
SBGs evolving out to Zm = 5.0, for different values of the evo¬ 
lution parameter. For each of the three values g = 0, 3, 4.5, indi¬ 
cated by solid, long-dashed, and dotted curves, the lighter curve rep¬ 
resents the unabsorbed bremsstrahlung emission plus absorbed GPL, 
and the heavy curve shows the sum of these two components plus 
Comptonized emission (which consists of an absorbed ai = 1.8 PL 
and unabsorbed a 2 = 2.3 PL). The local SBG X-ray luminosity den¬ 
sity is given in eq.(4), its cosmic evolution is given in eq.(5), and the 
template SBG spectrum is shown in Fig.2. The solid black curve on 
top represents the observed CXB [see eq.(ll)]. 

of the pure luminosity type (Lilly et al. 1995; Ellis et al.l996; 
Fried et al. 2001), or of the pure density type (Lin et al. 1999), 
is deduced also from optical surveys. 

The luminosity density of SBGs evolves mainly as result 
of their steeply varying density, and - to a smaller extent, as 
discussed below - the changes in the relative strength of the 
spectral components. Density evolution is represented by the 
factor (1 ± z)‘^, with q > 3 after the process of galaxy forma¬ 
tion began in ernest at some early redshift, Zm- Since we do 
not know when exactly galaxies began forming, the modelling 
the starburst phase of early galaxies necessitates viewing both 

(Franceschini et al. 2002) suggests properties fairly similar to those 
of SBG spectra: i.e., the ‘starburst signature’ discussed in section 2, 
plus sometimes a hard PL component (interpretetd as SB-obscured 
AGNs). ULIRGs are then implicitly included in the calculations re¬ 
ported in this paper. 
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q and Zm as essentially free parameters. ^ We factor out the 
dominant redshift dependence of the luminosity density: 

C{z) = C{Q) {I + zf , (5) 

with all the spectral dependence on z collected in the function 
this function will be specified below. The evolution of 
SBGs is then largely gauged by the index q an estimate for 
£(0) is given in eq.(4). 


5. SBG Emission as a Background 

The surface brightness due to a population of unresolved 
sources is given by (e.g., Boldt 1987) 

/ = ( 6 ) 


where C = J L dL is the comoving volume emissivity, 
and d£ is the line element. Since d£ is related to the look-back 
time r and the redshift zbyd£ = c/(l + z)dT, we obtain 



The ensuing spectral density is then given by 

^ ± [ ^dr, 

de 47r J deo ’ 


(7) 

( 8 ) 


where e is the energy measured in the observer frame and cq = 
e (1 + z) the energy measured in source rest frame. 

The look-back time r is related to the cosmic time t by 
T = to — t (with to the present age of the universe), and t is 
related to the redshift z by 


df 11 

dz 1 -I- z i7o 


-f 17^(1 + -f T (1 — ^o) (1 T z)"^ 


^ As an example, we derive an evolution function of the type in 
eq.(5) within a toy model of galaxy evolution. The comoving number 
density of galaxies undergoing a starbursting phase at any given epoch 
z, nsBG(z), is closely related to the total galaxy density at that time, 
no (z). If so, we expect the comoving SBG fraction to increase steeply 
with z. The interaction rate between galaxies is oc nG{z)^v (with 
constant interaction cross section), where naiz) is the galaxy density, 
and V = a/ (5 v)^ is the rms pairwise velocity dispersion of interacting 
galaxies. Because of the strong no oc (1 -I- z)® dependence we can 
ignore the relatively weak dependence on v. This implies that in a flat 
universe the volume density of SBGs at any given time is nsBG (z) = 
JrsBG(0)(l + z)"‘-®. 

® Evidence for evolution of the galaxy luminosity density has been 
found also in other bands. Based on z < 1.5, B- and f7k-selected 
galaxies, the deduced evolution rate has q — 2.4 ± 1.0 for an EdS 
cosmology and g = 1.7 ± 1.0 for a flat A model (Wilson et al. 2002). 
Earlier work, based on z < 1, /-band selected galaxies, claimed 4 
(Lilly et al. 1996). 




Fig. 4. The fraction of the CXB arising from SBGs as a function of 
energy, for various strengths of cosmological evolution. Models, pa¬ 
rameter values, and symbols are as in Fig.3. 


where fir, ^m, ^A, and flo are the current density parameters 
for radiation, matter. A, and their sum, respectively. In the mat¬ 
ter dominated era, 

J (l + z)-2-5_ _ EdS 

dz Ho \ [(1 -I- z)a/0.3 (1 -I- z)3 -I- 0.7] ^ flat A 

The SBG photon spectrum we have used in our calculations 
is given in eq.(l). Whereas the first two spectral components in 
eq.(l) have been adequately discussed in many previous works, 
the Compton CMB component needs to be further elaborated 
upon here because of its enhanced role in early SBGs. Most 
of the X-ray emission from a SBG is triggered by stellar activ¬ 
ity and involves the internal properties of the galaxy during a 
starburst phase of its evolution. As such, the emission is largely 
independent of redshift, if the main characteristics of this phase 
have not changed over cosmic time. The Compton component 
of the emission has evolved significantly, however, due to the 
much higher energy density of the CMB at z > 1. As we have 
noted in section 2, in a local SBG the FIR energy density is typ¬ 
ically much higher than that of the CMB, by as much as a factor 
of ~ 10 (Goldshmidt & Rephaeli 1995), so the Compton CMB 
component begins to dominate already at ^ 1. This, and the fast 
evolving SBG density, result in a very steep combined spectral 
and density evolution of the emission from SBGs at high X-ray 
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energies where Compton scattering is the main emission com¬ 
ponent in SBGs. Due to the steep (1 -|- z)^ dependence of the 
CMB energy density, z) has an explicit dependence 

on z in addition to the implicit dependence through x. 

With £(eoj z) = C{z) (j){eo), where £{z) is given in eq.(5) 
and (j>{€o) (X eg /(eo) is the SBG energy spectrum (source 
frame) normalized by the 2-10 keV flux, we have for the total 
intensity at an observed energy e due to a population of SBGs 
undergoing cosmic evolution according to eq.(5) out to a max¬ 
imal redshift Zm 


m 


c£{0) 

inHo 



(i+zY4>[t{i+z),z] ^ ^ 

--tlz 

(l+zy^\e{l+z),z] 

(r-|-2)^0.3(l+2)3-|-0.7 


EdS 

dz flat A 


( 10 ) 


The integrated intensity from SBGs was calculated at (the 
source frame) energies eg = 2-120 keV, for values of the evo¬ 
lutionary index q in the range 0.0 < g < 4.5, and adopting 
Zm = 5. In Fig.3 we show the integrated background light aris¬ 
ing from an evolving population of SBGs as a function of the 
(observed) energy, e. To determine the fractional contribution 
of SBGs to the CXB intensity in the 2-20 keV band, for the 
latter we use the analytic fit: 

/cxB(e) = 7.877 (11) 


with /cxB measured in keV/ (cm^ s keV sr) and e mea¬ 
sured in keV, given by Gruber et al. (1999) in the 3-60 keV 
band (based on HEAO-1/A2+A4 all-sky data: dominated by 
A2 proportional-counter data for e < 15 keV, and by A4 
scintillator-experiment data for e > 30 keV) and shown by 
Gendreau et al. (1995) and Chen et al. (1997) to hold down to 1 
keV (from ASCA data, based on a composite sample of the sky 
corresponding to less than a square degree and therefore lim¬ 
ited by by sky surface brightness fluctuations). ^ In Fig.4 we 
plot the SBG fractional contribution to the CXB as a function 
of (observed) energy. 

Looking at Fig.3 we see that for energies 5 keV the 
shape of the SBG background light gets progressively steeper 
than that of the CXB, for the range of values of q considered 
here. (The steep rise of the Comptonized CMB contribution 
moderates the decline somewhat.) From Fig.4 it is seen that for 

’ Although the spectral shape of the CXB is very well established, 
there is some question about the normalization, especially at 3 keV. 
Under the assumption that Gruber et al.’s (1999) composite multi¬ 
experiment fit [which in the 3-50 keV band is very close to the spec¬ 
trum implied by the HEAO-UAl data alone (Marshall et al. 1980; 
Boldt 1987)] is valid down to 1 keV, we note that the HEAO-1 spec¬ 
trum is uniformly lower, by ~ 11%, than the PL fit to the ASCA/SIS 
spectrum, 7cxB(e) = 8.6 keV/(cm^ s keV sr) (Gendreau et 

al. 1995). Actually, the discrepancy could be even worse, were the 
HEAO-1/A2 result to be renormalized downwards to account for the 
fact that it includes relatively bright unresolved sources that would 
have been resolved by ASCA. In addition, the ROSAT observed CXB 
(at £ > 1 keV), though limited by surface brightness fluctuations esti¬ 
mated at ^ 10% (Miyaji et al. 1998), suggests a CXB normalization 
^ 20% higher than that obtained from the HEAO-1/A2 result. Given 
this state of affairs, we assume that a smooth extrapolation of Gruber 
et al.’s (1999) HEAO-l-based formula down to 2 keV is safe enough 
for our current purposes. 


no cosmic evolution (q = 0) the contribution of SBGs to the 
CXB is < 1% at all energies; while for an evolution as steep 
as the maximum permitted by IR data (q ~ 3) the contribu¬ 
tion is ~ 5% at 5 keV, ~ 2% at 10 keV, ~ 1% at 15 keV, and 
vanishing at higher energies. This is so for the EdS model; the 
corresponding results are even lower for the flat A model. 

Also of interest is to determine the z dependence of the 
superposed SBG emission (at a given energy) in order to assess 
when the contribution is maximal. This can be readily found 
from the integrand of the z-integrated spectral intensity in eq. 
(10). Representing the main part of the spectral energy flux of 
a SBG as (()(eo) oc eo/(eo) cx with T ~ 2 in the 1- 

200 keV band (see Fig.2), we see that in the EdS model the 
integrand in eq. (10) is oc (1 -f z)‘^~^'^~^. Thus, in evolving 
models the relative contribution of SBGs increases with z if 
q > 3.5, and decreases if g < 3.5. While the critical value of g 
for which this transition occurs is different in the flat-A model, 
the behavior is qualitatively similar. Note that, as emphasized 
already, the contribution of the Compton component of the flux 
grows much more steeply due to the fact that the CMB energy 
density is oc (1 -I- z)^. 

6. Discussion 

The starting point in our treatment and previous similar works 
is the identification of SBGs based on their spectral proper¬ 
ties at 60/im, particularly the use of the same estimate for the 
local 60/im luminosity density, Ceo/im- However, most previ¬ 
ous studies were limited to the consideration of only the low- 
energy (<2 keV) contribution of SBGs to the CXB. Weedman 
(1987), Griffiths & Padovani (1990), and Ricker & Meszaros 
(1993), having based their estimates on Einstein data, evalu¬ 
ated the SBG contribution at 2 keV. More recently, Moran et 
al. (1999) estimated the SBG contribution to the 5 keV CXB 
from ASCA data, based on the assumption that the 5 GHz ra¬ 
dio emission from galaxies lacking a radio-loud AGN can serve 
as a gauge of star-formation activity, and that the value mea¬ 
sured in NGC 3256 for the 5 keV to 5 GHz flux ratio is univer¬ 
sal. They used the 5 GHz number-count-flux relation, logA^- 
logS'sGHz, in the sub-mJy regime, and estimated that SBGs 
contribute ~ 12% — 28% of the 5 keV CXB. A less limited as¬ 
sessment of the SBG contribution was implemented (Rephaeli 
et al. 1991, 1995) by considering the SBG emission in the 2-30 
keV band. However, due to lack of observational data, our pre¬ 
vious works were based on a very preliminary statistical anal¬ 
ysis of stacked (low-exposure HEAO-1 A2h-A 4) data. Here, on 
the other hand, we have estimated the 2-10 ke V luminosity den¬ 
sity using the detailed 60/im and 2-10 keV luminosities mea¬ 
sured for essentially the same sample of local SBGs whose X- 
ray spectra have been adequately measured (though not yet at 
the desired level of precision). 

Inclusion of nonthermal emission due to Compton scatter¬ 
ing of relativistic electrons by the FIR and CMB fields is an im¬ 
portant new ingredient in the modeling of the high energy emis¬ 
sion in an evolving population of SBGs. Although the level of 
this emission is still unknown, our explicit accounting for this 
emission and its expected rapid rise with redshift (which makes 
its relative contribution very significant at z > 1) provide fur- 
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ther motivation for a more detailed description then given here. 
[The occurrence of low-luminosity AGNs in the nuclei of SBGs 
is a controversial issue (e.g., Veilleux 2001), so we chose not 
to include such a component in the current treatment.] 

Our main result that the evolving population of SBGs con¬ 
tributes appreciably to the CXB in the 2-10 keV band is quite 
consistent with our previous estimates (Rephaeli et al. 1991, 
1995), and a similar result by Griffiths & Padovani (1990). It 
was argued in the latter paper that since the spectrum of X- 
ray binaries is quite hard and comparable in shape to that of 
the residual CXB in the 2-30 keV range, evolving SBGs in the 
redshift range 0.5^ 1 could make an important contribu¬ 

tion to the 2-10 keV CXB. In the use of a quantitative spec¬ 
tral form and evolutionary function, our treatment is similar in 
spirit to that of Ricker & Meszaros (1993, hereafter RM93). 
The main improvement with respect to the work of RM93 is 
the use of spectral results from individual SBGs observed with 
ASCA and BeppoSAX and the inclusion of the two Compton 
components, in contrast to their use of essentially an assumed 
theoretical model. [We have included also the observed warm 
(~ 0.8 keV) component, but its contribution to the emission at 
e > 2 keV is relatively minor.] 

In order to check the stability of our results vs. the char¬ 
acteristics of the FIR-selected SBG population, we have re¬ 
peated the calculation using the SBG population parameters of 
Franceschini et al. (2001), who fitted the 12/im number counts 
invoking the presence of a starbursting population whose vol¬ 
ume emissivity increases steeply out to z = 0.8, and remains 
constant at higher z [in a flat A (fim = 0.3, Ha = 0.7) 
cosmological model]. The contribution of Franceschini et al.’s 
SBG population * turns out to be only marginally higher than 
that computed for the SBG population of Pearson & Rowan- 
Robinson (1996) used in our main calculation and described by 
eqs.(2)-(5) (with q = 3). On the whole, then, the results of this 
variant calculation are qualitatively consistent with our main 
results (see Fig.5, where the ’wind-tbinaries’ spectral model is 
used). 

Recent results from the Chandra deep surveys provide di¬ 
rect estimates of the SBG contribution to the CXB. Using the 1 
Ms Chandra Deep Field North (Brandt et al. 2001a) and 15 p,m 
ISOCAM Hubble Deep Field North (Aussel et al. 1999) sur¬ 
veys, Alexander et al. (2002) found a tight correlation between 
the population of strongly evolving SBGs at z ~ 1 discovered 
in faint 15p,m ISOCAM surveys (e.g., Aussel et al. 1999; Elbaz 
et al. 1999, 2002) - whose space density is an order of magni¬ 
tude higher than that of local SBGs - and the apparently normal 
galaxies detected at faint fluxes in X-ray surveys (Giacconi et 

* The local 2-10 luminosity density of the starbursting population 
of Franceschini et al. (2001) is £(0) = 3.15 x 10®® erg s“® Mpc 

This value is deduced from the 12/im local luminosity func¬ 
tion of this population (see their Fig.9), which can be described as 
<(.(L)dlogL = <)-* L^““[l-fL/(6L*)]"'’dlogL with = 5.5 X 
10“®, L* = 5 X 10®Lq, a = 1.15, and 6 = 3.1, and the average 2-10 
keV to 12/im luminosity ratio of (1.26 ±0.35) x 10“®, derived for the 
same sample of local SBGs mentioned in footnote 3 (the 12/im data 
are from Shapley et al. 2001). The evolution is £(z) = £(0)(l±z)^'® 

for Z < Zbreak and C{z) = £(0)(1 ± Zbreak)’’ ® for Z > Zbreak, with 
Zbreak ~ 0.8. 



e [keV] e [keV] 

Fig. 5. The spectrum of the integrated background light (left), and the 
corresponding contribution to the CXB (right), arising from the SBG 
population of Franceschini et al. (2001), assuming a ’starburst only’ 
model for the spectrum. See the text for more details. 

al. 2001; Hornschemeier et al. 2001; Brandt et al. 2001b). The 
characteristics (Lfir, L^, 2-10 keV spectral slope) of the pop¬ 
ulation of the X-ray detected galaxies, whose redshifts are in 
the range 0.4 < z < 1.3 (based mainly on Keck data; see 
Alexander et al. 2002 and references therein), are consistent 
with those expected for SBGs resembling local objects such 
as M82 and NGC 3256. In particular, their stacked average 2- 
10 keV spectral slope, F ~ 2, is steep enough to suggest a 
low fraction of obscured AGN activity within the population. 
The contribution of this population to the 0.5-8 keV CXB, es¬ 
timated to be ~ 2% (Alexander et al. 2002), marginally agrees 
(on the low side) with our result for q = 3 (see Fig.4). A con¬ 
tribution to the CXB will also come from the 1 Ms Chandra 
Deep Field North sources identified as Lyman-break galaxies 
at z ~ 2 — 4: these star-forming sources have rest-frame 2-8 
keV luminosities and X-ray/i?-band luminosity ratios that are 
comparable to those of local SBGs (Brandt et al. 2001c). Based 
on the above, there is no evidence for a substantial change in 
the ~ 2 — 10 keV luminosity and spectral slope of these X- 
ray-selected, star-forming galaxies. On the other hand, their IR 
(15/iTO ISOCAM) counterparts do show a strong density evo¬ 
lution with redshift (e.g., Aussel et al. 1999; Elbaz et al. 1999, 
2002). All this clearly supports our adopted scheme of cosmic 
evolution, which is essentially of the density type (see sect.4). 

Einally, it is commonly thought that AGNs are the prime 
sources of the CXB (e.g.: Comastri et al. 1995; Setti & Woltjer 
1989; Boldt & Leiter 1981 and Leiter & Boldt 1982), and that 
additional contributions from other sources are typically esti¬ 
mated to be ^ 20% in the 2-10 keV range (Giacconi et al. 2001; 
Tozzi et al. 2001; Mushotzky et al. 2000). Based on our results, 
SBGs may possibly be the second most important contributors 
to the CXB. (Clusters of galaxies also contribute a few percent: 
e.g., Piccinotti et al. 1982.) The importance of determining the 
level of the SBG contribution to the 2-10 keV CXB lies not 
only in the implications for the evolution of SBGs themselves, 
but also in the ramifications for galaxy evolution in general, and 
as an additional input on the evolution of the AGN phenomenon 
in particular (e.g.: Eabian et al. 1998; Gilli et al. 1999). 
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